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ABSTRACT 

Star formation in relic HII regions of the first stars is investigated using 
magneto- hydrodynamical simulations with a nested grid method that covers ~ 10 
orders of magnitude in spatial scale and ~ 20 orders of magnitude in density 
contrast. Due to larger fraction of H 2 and HD molecules, its prestellar thermal 
evolution is considerably different from that in the first star formation. Reflecting 
the difference, two hydrostatic cores appear in a nested manner: a protostar is 
enclosed by a transient hydrostatic core, which appears during the prestellar col- 
lapse. If the initial natal core rotates fast at a rate with rotational to gravitational 
energy ratio (3q > 0.01 —0.1, the transient hydrostatic core fragments to ~ 10 M® 
sub-cores at density ~ 10 9 cm -3 . With smaller rotation energy, fragmentation 
occurs at higher density while a single protostar forms without fragmentation if 
rotation is extremely slow with (3q < 10~ 6 — 10 -5 . If magnetic field is present, 
these threshold values of /3 is boosted owing to angular momentum transport 
by the magnetic breaking. Magnetic field also drives the protostellar outflows. 
With strong magnetic field, two distinct outflows are observed: The slower one 
emanates from the transient hydrostatic core, while the faster one from the pro- 
tostar. These flows may affect the final stellar mass by ejecting some of masses 
in the initial core, and also may play some role in driving and maintenance of 
interstellar turbulence in young galaxies. 

Subject headings: cosmology: theory — galaxies: formation — stars: formation 

1. Introduction 

Theoretical studies in the last decade revealed that first stars in the universe were born 
in redshift z ~ 20 — 30 as very massive stars with > 1OOM (e.g., Ciardi & Ferrara 2005). 
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It is also suggested that a fraction of them are formed as binaries or multiples as long as 
their natal cores somewhat rotate (Machida et al. 2008a; Turk, Abel, & O'Shea 2009). To 
drive protostellar outflows, ubiquitously observed in the solar neighborhood and affecting the 
protostellar accretion, from the first stars, however, requires a magnetic field of 10~ 8 — 10 _9 G 
at 10 3 cm -3 (Machida et al. 2006a), higher than the level generated by mechanisms currently 
proposed, e.g., Biermann batteries at the cosmic recombination or halo virialization (Ichiki 
et al.; Xu et al. 2008) by 2-3 orders of magnitudes. Hence the outflows/jets are hard to 
emerge at least in the early phase of protostellar accretion. 

Even in the primordial gas, prestellar thermal evolution, which affects fragmentation 
and thus characteristic stellar mass, can differ in a cloud disturbed by a former generation 
of stars. The stars born in such an environment begins to be recognized as a different class 
of objects and called Population III. 2 stars, while the genuine first stars are Population III.l 
stars (O'Shea et al. 2008). Various feedbacks can affects cradles of subsequent generation 
of stars, and thus Pop III. 2 in principle includes distinct sub-classes. The most studied 
example, though, is those formed in relic HII regions of the first stars. 

A large HII region forms around a massive first star. Pop III stars with mass > 100M Q 
end their lives as supernovae only in a limited mass range (140 — 26OM ; Heger & Woosley 
2002; Umeda & Nomoto 2002), and thus the most of them would collapse directly as black 
holes rather quiescently. The surrounding gas is left ionized and starts recombining gradually, 
and another episode of star formation commences after ~ 100 Myrs. In such relic HII regions, 
abundant electron, which is the catalyst for H 2 formation, allows the rapid replenishment 
of H2 to the fraction ~ 10~ 3 , about an order of magnitude higher than in the Pop III.l 
case. The cooler environment (< 150K) due to the higher H 2 cooling is preferable for HD 
formation, and its cooling further lowers the temperature down to a few tens of Kelvins, 
close to the CMB temperature (Nagakura & Omukai 2005). The fragmentation mass-scale 
due to the HD cooling is about ~ 40M Q , and most of it is expected to be accreted by the 
formed star eventually (Yoshida, Omukai, & Hernquist 2007; McGreer & Bryan 2008). Not 
only in the relic HII regions, but the HD cooling can be important also in a medium ionized 
at the virialization of a large halo or even in a small halo if dynamical heating during mass 
assembly is weak (Ripamonti 2007). This "HD mode" of primordial star formation may be 
a more common mode than the Pop III.l one in forming galaxies. 

In addition, if all the Pop III.l stars were indeed very massive with a few hundreds of 
solar mass, they might have been very inefficient cauldrons for metals owing to the limited 
mass range for supernovae. With smaller mass and being able to explode as Type II super- 
novae, Pop III. 2 stars could be more responsible for the metals in the early universe. To 
predict its mass range more accurately, we need to know the conversion efficiency from a 



- 3- 



dense core to stars, which is determined by binarity/multiplicity, feedback from protostellar 
outflow etc. As an approach to this problem, here, we study the prestellar collapse of rotat- 
ing dense cores that is produced in an initially ionized gas, putting main emphasis on the 
condition of binary formation and launching of outflows. 

This paper is organized as follows. The numerical method is briefly described in §2, 
and the results are presented in §3. Finally, we discuss its significance in the cosmological 
context and caveats in §4. 

2. Numerical Method 

First we calculate the temperature evolution by using a free-falling one- zone model of 
Omukai et al. (2005), starting from number density 1cm -3 , temperature 7000K and the 
initial ionization degree 0.1. The CMB radiation at 27.3K is imposed. With given temper- 
ature evolution as a function of density, we performed magneto-hydrodynamics (MHD) by 
a nested-grid method to gain high resolution near the center. As the initial condition, we 
take the critical Bonner-Ebert sphere with central density 10 4 cm -3 but with a factor 1.68 
of density enhancement to induce the gravitational instability. In addition, m = 2 mode of 
density perturbations are imposed with 10% to break axi-symmetry and cause fragmenta- 
tion. The initial angular momentum and magnetic field are parallel to the z-axis, and are 
parameterized by the ratio of rotation to gravitational energy, /3 , and the ratio of magnetic 
to gravitational energy 70, respectively. These quantities are related to the initial rotation 
rate of the core f2 (s _1 ) and magnetic flux density B (G) as (3 — 0.1[Q /3 x 10" 14 (s -1 )] 2 
and 70 = 1.5 x 10~ 5 [.Bo/iO~ 7 (G)] 2 , respectively. The magnetic field is assumed to be frozen 
to the gas. For details of the MHD and initial conditions, see Machida et al. (2005). The 
initial parameters of calculated runs and summary of the result are presented in Table 1. 

3. Results 

The adopted temperature evolution in the HD mode is shown in Figure 1. Also shown 
is that in the Pop III. 1 case ("H 2 mode"), where the collapse begins from an almost neural 
condition. The temperature in the HD mode is lower than that in the H2 mode for 10 2_7 cm -3 , 
and reaches several 10K owing to the HD cooling. The HD-cooling rate saturates at ~ 
10 5 cm~ 3 by collisional de-excitation, and the temperature begins to increase with density. 
The fragmentation scale ~ 40M Q is set around there (Yoshida et al. 2007). At 10 8 cm~ 3 , the 
H 2 begins to form via the three-body reaction, whose heating makes the temperature jump 
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up abruptly. The simultaneous increase of pressure slows down the gravitational collapse 
near the center of the dense core, and transient hydrostatic core forms just for a brief period. 
The core begins to collapse dynamically again until the permanent hydrostatic core, i.e., the 
protostar, forms at > 10 21 cm~ 3 unless the centrifugal supports prevents the collapse before 
that. Similar transient hydrostatic objects, the so-called first cores, form in metal-enriched 
cores when their central parts become optically thick to the dust absorption. 

In metal-enriched cases, the first-core formation is an important epoch for fragmentation 
of rotating cores (Machida et al. 2009). In the isothermal phase, i.e., before the first-core 
formation, even though the core becomes flattened by rotation, its collapse does not signifi- 
cantly slow down until the adiabatic temperature-increase leading to the first-core formation. 
Shortly after the first-core formation, if the core is flat enough, the slow contraction allows 
perturbations to grow, thereby causing fragmentation into a binary/multiple system. On the 
other hand, in the H 2 mode of primordial star formation (i.e., PopIII.l formation), where 
the temperature continues to increase gradually without transient-core formation, no char- 
acteristic epoch for binary/multiple formation exists. In this case, when a thin disk forms 
by rotation, it fragments after further contraction of a few orders of magnitude in density. 

The final states of the calculations are presented for different combinations of rotation 
(/?o) and magnetic field strength (70) in Figure 2. The cores with initially large rotation 
(A) ^ 0.01 — 0.1 in the 70 = cases) fragment shortly after the transient core formation 
and form wide proto-binaries with initial masses ~ 1 — 1OM and separations ~ 10 4 AU 
(cases a in Figure 2). Those with smaller rotation (10~ 5 < /3q < 10 -2 in the 70 = cases) 
fragment at higher density before the protostar formation unless the rotation is extremely 
small (Po ^ 10~ 6 — 10~ 5 in the 70 = cases). Such high-density (10 17 ~ 20 cm~ 3 ) fragmentations 
produce close proto-binaries with initial masses ~ 1O~ 2 M and separations ~ 10~ 2 AU (cases 
b in Figure 2). Note that some models classified in case b show merger of fragments; i.e., 
models (fi , 70) = (10~ 2 , 1.5 x 10~ 5 ) and (10~ 4 , 1.5 x 10 -5 ). We expect that the merged cores 
finally evolve in a single massive star, since spiral structures in them effectively transfer the 
angular momenta outward, which prevents further fragmentation. Those with even smaller 
rotation collapse all the way without fragmentation and form single protostars (cases c in 
Fig. 2). While the collapse is approximately spherical, the rotation parameter (3 at the 
center of the core increases as oc n 1 / 3 , where n is the number density, from the initial value 
Po- When (3 reaches 0.04-0.1, the radial contraction is retarded by the centrifugal force. The 
core becomes a thin disk after more contraction by a few orders of magnitude in density 
(e.g., Machida et al. 2008). The threshold (3 values (i.e., 0.01 - 0.1 and 10~ 6 - 10~ 5 ) for 
binary formation correspond to the conditions of thin-disk formation before formation of the 
transient core (at ~ 10 8 cm~ 3 ) and the protostar (at ~ 10 21 cm~ 3 ), respectively. 
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The magnetic field reduces the angular momentum of the cores through magnetic break- 
ing, hence boosting the threshold rotation rate for fragmentation. For example, with initial 
magnetic field 70 = 1.5 x 10~ 3 (1.5 x 10 _1 ), the threshold f3$ for fragmentation increases 
above 10~ 4 (10~ 2 , respectively) from ~ 10~ 5 in the 70 = case. 

Although the initial mass of fragments is very small, they are embedded in a dense 
core and subsequently grow by accretion until the majority of the core material (< 1OM ) 
is accreted by the binary/multiple system (Yoshida et al. 2008). However, if launched, 
protostellar outflow may reduce the conversion efficiency from the core to the system by 
evacuating some of envelope materials. 

When the initially ordered magnetic field is highly twisted by rotation, outflows are 
launched either by magnetic centrifugal force or by pinching by magnetic pressure. During 
dynamical collapse, the cores do not rotate so many times as to twist the magnetic field 
enough. After hydrostatic core formation, the magnetic field is amplified and distorted 
owing to the slow contraction, thereby enabling the launching of slow outflows. In the case 
of present-day star formation, two kinds of outflows appear, corresponding to two epochs 
of hydrostatic core formation, i.e., first core and protostar (Machida et al. 2008b). In our 
case, owing to the presence of the transient hydrostatic core, similar two flows are launched, 
i.e., the slower one (2-3km/s) at transient core formation, and the faster one (20-40km/s) 
just at the protostar formation. The speeds of the flows increase with the mass of driving 
objects, i.e., the transient hydrostatic core and protostar. According to our experiments, 
the minimum initial magnetic field needed for driving the fast (slow) flow was about 10 _8 G 
(10~ 7 G, respectively). On the other hand, in the Pop III. 1 case, only fast flow is launched 
as no transient hydrostatic stage exists (Machida et al. 2006a). 



4. Discussion 

Currently, we have no information about rotation rates of Pop III. 2 dense cores. If we 
assume that they are similar to the values in the Milky Way (/3q ~ 10 -4 — 0.07 with the 
median of 0.02; Caselli et al. 2002), fragmentations typically occur at far higher densities 
than the transient core formation (i.e., case b in Figure 2), and produce close binary/multiple 
systems. Since the mass scale of the natal cores is 30 — 4OM , set by fragmentation induced 
by the HD cooling (Yoshida et al. 2007), the final product would be a close binary of 
10 - 2OM stars. 

For the Pop III.l protostars, only the primordial magnetic field of the cosmological origin 
is available. On the other hand, the Pop III. 2 stars are able to use additional magnetic fields, 
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which might have been generated and amplified by activities of previous generations of stars, 
including both the Pop III. 1 and III. 2. Therefore, the Pop III. 2 protostars could launch the 
outflow even though an order-of-magnitude larger magnetic field is required for the fast-flow 
launching than in the Pop III. 1 case. To know whether the protostellar jet is launched, we 
need to wait future numerical simulations to narrow down the possible values for magnetic 
field in such an environment. If indeed launched, the jet expels some materials in the natal 
core, and hence reducing the final mass in the binary system. Matzner & McKee (2000) 
estimated the efficiency from the core to the system at 25 — 70% after the jet feedback in the 
present-day star formation by way of hydrodynamical simulation. If this efficiency is similar 
in our case, the final system mass would be 2OM out of a 4OM core. To know the final 
mass of stars, we need to follow evolution in the accretion phase, which is left for future 
studies. 

As mentioned in §2, we used the temperature evolution as a function of density, which 
is calculated by the one-zone model of Omukai et al. (2005), where spherical clouds with 
the Jeans lengths are assumed to collapse at about the free-fall rate. In reality, thermal 
evolution depends on the cloud dynamics, and vice versa. For example, when the cloud 
shape deviates from the sphere, the optical depths and thus the cooling rates can differ from 
the assumed values. Also, slower collapse leads to lower compressional heating. Moreover, 
our method cannot capture the shock heating properly. This could be a significant drawback 
in simulating more turbulent situations (see Clark et al. 2008 for more discussions). 

The first galaxy, harboring the first star clusters, is found to be turbulent owing to 
the cold accretion flow at its formation by numerical experiments (Wise & Abel 2007; Greif 
et al. 2008). Here, however, we have not use turbulent initial conditions for simplicity. 
If the turbulent energy is comparable to the gravitational one, the cloud evolution can be 
significantly affected in particular in early phases leading to formation of dense cores (e.g., 
Clark et al. 2008). Rotational energy of the dense cores would be set by the turbulence and 
tend to be higher in the more turbulent environment. In this study, we have regarded the 
degree of rotation as a free parameter and followed the evolution only after the formation of 
the dense cores. 

The protostellar outflows by Pop III. 2 stars could contribute to drive and maintain the 
turbulence in the first galaxies. However, this phase might be short-lived as the first galaxies 
will be metal-enriched before long. In addition, such metal-free first galaxies might be rare as 
the previous generations of stars have already polluted with metals. Therefore, protostellar 
outflows by the Pop II stars would be more important in this context and awaits future 
studies (e.g., Machida et al. 2009). 
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Table 1: Models and Results 
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a Number density at fragmentation epoch 

a Emergence of Slow, Fast and Both (slow and fast) outflows. 
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n (cm- 3 ) 

Fig. 1. — Thermal evolution of a metal- free prestellar core from an initially ionized state 
as function of the number density (thick solid; indicated as "Pop III.2"). The temperature 
divided by mean molecular weight fi is shown. Also shown for comparison is that for the Pop 
III. 1 case (dashed). The lines of constant Jeans mass are presented by the thin dot-dashed 
line. Some characteristic epochs are indicated by arrows with captions. 
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Fig. 2. — The final states for runs with different rotation ((3q) and magnetic field strength (70)- 
The results are classified according to fragmentation properties; (a) fragmentation shortly after the 
transient core formation, (b) fragmentation at higher density, and (c) no fragmentation, which are 
indicated by the background colors. Note the larger binary separations (~ 10 AU) in the cases (a) 
than those in (b; ~ 10 -2 AU). In each panel, the structure of high-density region (red isodensity 
surface) and magnetic field lines (black and white streamlines), and outflow shape (the region of 
radial velocity v r > 0; transparent orange surface) are plotted in three dimensions, while the density 
contour (color and contour lines) and velocity vectors (arrows) are projected on each wall surface. 
The grid level I and grid scale are shown at the upper-left corners and bottom, respectively, of each 
panel. 



